DNA nucleotide-specific modulation of fiA transverse edge currents through a 
metallic graphene nanoribbon with a nanopore 
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We propose two-terminal devices for DNA sequencing which consist of a metallic graphene 
nanoribbon with zigzag edges (ZGNR) and a nanopore in its interior through which the DNA 
molecule is translocated. Using the nonequilibrium Green functions combined with density func- 
tional theory, we demonstrate that each of the four DNA nucleotides inserted into the nanopore, 
whose edge carbon atoms are passivated by either hydrogen or nitrogen, will lead to a unique change 
in the device conductance. Unlike other recent biosensors based on transverse electronic trans- 
port through DNA nucleotides, which utilize small (of the order of pA) tunneling current across a 
nanogap or a nanopore yielding a poor signal-to- noise ratio, our device concept relies on the fact 
that in ZGNRs local current density is peaked around the edges so that drilling a nanopore away 
from the edges will not diminish the conductance. Inserting a DNA nucleotide into the nanopore 
affects the charge density in the surrounding area, thereby modulating edge conduction currents 
whose magnitude is of the order of /jlA at bias voltage c± 0.1 V. The proposed biosensor is not lim- 
ited to ZGNRs and it could be realized with other nanowires supporting transverse edge currents, 
such as chiral GNRs or wires made of two-dimensional topological insulators. 

PACS numbers: 87.14.G-, 73.63.Rt, 72.80.Vp 



The successful realization of fast and low-cost methods 
for reading the sequence of nucleotides in DNA is envis- 
aged to lead to personalized medicine and applications in 
various subfields of genetics. The solid-state nanopores 
represent one of the pillars of the so-called third gen- 
eration sequencing [Tj. The key issues in this approach 
revolve around how to slow down the translocation speed 
of DNA and how to achieve single-nucleotide resolution. 

The very recent experiments [2]-[4] on DNA transloca- 
tion through graphene nanopores have introduced a new 
contender into this arena. Graphene — the recently dis- 
covered [5] two-dimensional allotrope of carbon whose 
atoms are densely packed into a honeycomb lattice — 
brings its unique electronic and mechanical properties 
into the search for an optimal nanoelectronic biosensor. 
Since single layer graphene is only one-atom-thick, the 
entire thickness of the nanopore through which DNA is 
threaded is comparable to the dimensions of DNA nu- 
cleotides. Therefore, there is only one recognition point 
rather than multiple contacts with DNA in the nanopore. 

However, the recent experiments [2-4 on nanopores 
within single or multilayer large- area graphene, which 
have measured fluctuations in the vertical ionic cur- 
rent flow due to DNA translocation through the pore, 
have not reached sufficient resolution to detect and iden- 
tify individual nucleotides. An alternative scheme is 
to adapt the transverse current approach to graphene- 
based biosensors [6H9]. The past several years have 
seen a number of theoretical proposals [lOj [IT] and ex- 
periments P2J [13] on nanogaps between two metallic 
electrodes where the longitudinally translocated DNA 
through the gap modulates the transverse tunneling cur- 
rent. Also, the recent first-principles simulations have 




FIG. 1: Schematic view of the proposed two-terminal device 
where transverse conduction current flows around the zigzag 
edges of a metallic graphene nanoribbon with a nanopore, 
while DNA molecule is translocated through the pore to 
induce nucleotide-specific-modulation of such edge currents. 
The device central region, which is simulated via NEGF- 
DFT formalism, consists of 14-ZGNR (composed of 14-zigzag 
chains which determine its width ~ 3.1 nm) and a nanopore of 
^1.2 nm diameter. The edge carbon atoms of the nanopore 
are passivated by either hydrogen or nitrogen, while edge 
atoms of ZGNR itself are passivated by hydrogen. The total 
number of simulated atoms (C-blue, H-yellow, N-green, O- 
red, P-orange) in the central region, including the nucleotide 
within the nanopore, is around 700. 



analyzed such modulation of the tunneling current for 
a nanogap [7J [8] between metallic GNRs with zigzag 
edges (ZGNR) or a nanopore [9 within semiconducting 
graphene nanoribbons with armchair edges (AGNR). 

However, the tunneling-current based graphene biosen- 
sors will face the same challenges [IT] encountered by cur- 
rent experimental efforts to utilize transverse tunneling 
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FIG. 2: (a) The room-temperature conductance of the two- 
terminal 14-ZGNRs with ^1.2 nm diameter nanopore whose 
edge carbon atoms are passivated by either hydrogen (H-pore) 
or nitrogen (N-pore). (b) The room-temperature conduc- 
tance of the same device as in panel (a) when one of the 
four DNA nucleotides (A-adenine, C-cytosine, G-guanine, T- 
thymine) is inserted into the center of the nanopore within 
the 2/z-plane [Fig. |3je)]. These conductances are computed 
via first-principles quantum transport simulations where both 
panels compare results obtained using two different NEGF- 
DFT codes— our home-grown MT-NEGF-DFT [H [15] and 
commercial ATK 1161. 



current across a gap between two gold electrodes [T2J [13] , 
such as poor signal-to-noise ratio at small bias volt- 
ages due to the fact that molecular eigenlevels are typ- 
ically far away from the Fermi energy of the electrodes. 
In this case, the tunneling is off-resonant and currents 
are of the order of pA at typically applied bias voltage 
— 0.5 V 112] [13]. Such small off-resonant tunneling cur- 
rents are highly dependent on dimcult-to-control relative 
geometry between the molecule and electrodes, so that 
recent experiments have measured broad current distri- 
butions corresponding to each nucleotide in the case of 
bare gold electrodes [12] and somewhat narrower but still 
overlapping distributions [13] for functionalized gold elec- 
trodes. Similarly, first-principles simulations of tunneling 
through the nanogap hosting a DNA nucleotide between 
two metallic GNRs have revealed current variation over 
several orders of magnitude (e.g., 10 -2 -10 -10 nA at bias 
voltage 1 V [7]) when changing the position and orienta- 
tion of nucleotides within the gap. 

The theoretical proposals to increase transverse cur- 
rent across the nanogap, as in the case of carbon nan- 
otube electrodes terminated with nitrogen where in- 
troduction of states closer to the Fermi level enables 
quasiresonant tunneling, offer only slight improvement — 
nA current at ~ 0.4 V bias voltage [11 . Applying higher 
bias voltage to increase the current signal is detrimental 
since it can lead to attraction of the negatively charged 
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FIG. 3: (a) The variation of the room-temperature conduc- 
tance of 14-ZGNR with N-pore due to the rotation of A, C, G, 
T nucleotides within the nanopore. The shaded vertical rect- 
angles mark the regions of overlap between the conductance 
intervals associated with different DNA nucleotides, (b)-(e) 
The specific positions of a nucleotide (guanine in the exam- 
ple) within the N-pore that define the conductance intervals 
shown in panel (a). The conductances in panel (a) were com- 
puted using our home-grown MT-NEGF-DFT code [T4l fT5] . 



DNA backbone toward one of the electrodes thereby im- 
peding the translocation. 

Here we propose a novel device concept that could re- 
solve these issues by abandoning the usage of small tun- 
neling current altogether. Its operation crucially relies 
on the existence of nanowires in which the spatial cur- 
rent profile [17] is confined around their transverse edges, 
so that drilling a nanopore should not change signifi- 
cantly their conductance which is of the order of few 
conductance quanta 2e 2 /h. When one of the four DNA 
nucleotides — adenine (A), cytosine (C), guanine (G), or 
thiamine (T) — is inserted into the nanopore in the course 
of DNA translocation, it will affect the charge density 
around the pore via its electrostatic potential thereby 
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modulating edge conduction currents that are several or- 
ders of magnitude larger than tunneling currents across 
nanogaps [6H3 HOI H2 H3] or nanopores [9] . 

The candidate nanowires supporting edge currents can 
be found among GNRs with zigzag edges or the very 
recently fabricated [18] chiral GNRs, as well as among 
two-dimensional topological insulators (2D TI) [19]. In 
the case of zigzag or chiral GNRs, spatial profile of local 
currents carried by electrons around the charge neutral 
point (CNP) shows large magnitude around the edge [20] 
and a tiny current flowing through their interior. In 2D 
TI nanowires, similar situation will appear if the wire is 
narrow enough so that helical edge states overlap slightly 
and edge currents can be modulated. Otherwise, in suffi- 
ciently wide 2D TI wires current is strictly confined to the 
edges and cannot be affected by time-reversal-preserving 
impurities, vacancies or modulation of charge density be- 
cause of the fact that helical edge states guide electrons 
of opposite spin in opposite directions to prevent their 
backscattering [T9] . 

The recent proliferation of nanofabrication tech- 
niques [18] [2TJ [22] for GNRs with ultrasmooth edges are 
making them widely available, and their exposed surface 
allows for an easy integration into biosensors. Therefore, 
in the device depicted in Fig. [T] we choose GNR with 
zigzag edges [22 . The device corroborates the general 
modulation-of-edge-currents concept discussed above, as 
demonstrated by our central result in Fig. [2] obtained 
via first-principles quantum transport simulations using 
two completely different [T4HT6] computational imple- 
mentations of the nonequilibrium Green function cou- 
pled to density functional theory (NEGF-DFT) formal- 
ism [23] [24]. 

Figure [2] shows how each DNA nucleotide inserted in 
the center of the nanopore within the ?/z-plane [Fig.[3je)] 
will change the device room-temperature conductance 
by a specific amount. When spatial orientation of nu- 
cleotides with respect to the pore is changed as in 
Fig. |3jb)-(d), the conductance will vary within the in- 
tervals shown in Fig. |3ja). The one-atom-thick GNRs 
make it possible to evade situations where several nu- 
cleotides inside the nanopore affect transverse conduc- 
tion simultaneously, as could be the case when mem- 
brane carrying the nanopore is insufficiently thin [25] . 
The DNA nucleotide-specific modulation of current / is 
achieved while remaining in the linear-response regime, 
where / = GV is of the order of fiA at bias voltage 
~ 0.1 V. Such sizable operating current is expected to 
be much larger than electrical noise due to thermal fluc- 
tuations of the DNA structure. 

In the NEGF-DFT formalism [23] [24], the Hamil- 
tonian is not known in advance and has to be com- 
puted by finding converged spatial profile of charge 
via the self-consistent DFT loop for the density ma- 
trix p = ^ J dEG < (E) whose diagonal elements give 
charge density [24] . The NEGF formalism for steady- 



state transport operates with two central quantities, re- 
tarded G(E) and lesser Green functions G < (E), which 
describe the density of available quantum states and how 
electrons occupy those states, respectively. In the co- 
herent transport regime (i.e., in the absence of electron- 
phonon or electron-electron dephasing processes), only 
the retarded Green function is required to post-process 
the result of the DFT loop by expressing the zero-bias 
electron transmission function between the left (L) and 
the right (R) electrodes as: 

T(E) = Tr {T R (E)G(E)T L (E)G\E)} . (1) 

The matrices T LiR (E) = i[E LiR (E) - ^ LjR (E)] account 
for the level broadening due to the coupling to the 
electrodes, where T>L yR (E) are the self-energies intro- 
duced by the ZGNR electrodes [24 . The retarded 
Green function matrix of the central region is given by 
G = [ES - H - £ L - EiJ -1 , where in the local orbital 
basis Hamiltonian matrix H is composed of elements 
Hij = (^il^Ksl^j) and Hks is the effective Kohn-Sham 
Hamiltonian obtained from the DFT self-consistent loop. 
The overlap matrix S has elements Sij = (<pi\<pj). 

The conductance at finite temperature T is obtained 
from the transmission function T(E) using the standard 
Landauer formula for two-terminal devices 

+oo 

G= J dET(E)(-^j, (2) 

— oo 

where f(E) = {1 + exp[(£ - ^/ksT}}- 1 is the Fermi 
function of the macroscopic reservoirs into which semi- 
infinite ideal leads terminate. The electrochemical po- 
tential \i is the same for both reservoirs at vanishingly 
small bias voltage. 

The retarded Green function G is computed for the 
central region finite- ZGNR+nanopore of the biosensor 
shown in Fig. [I] consisting of around 700 atoms. This 
central region is attached to two semi-infinite ZGNRs 
electrodes of the same width. Whereas graphene is me- 
chanically strong, it can be used as both the membrane 
material carrying a nanopore and the electrode material. 
In real devices, ZGNR electrodes will eventually need to 
be connected to metallic electrodes attached to an exter- 
nal battery. However, the fact that GNRs used in ex- 
periments are typically rather long and screening takes 
place over a distance much shorter [24] than the central 
region justifies the usage of semi-infinite ZGNRs as two 
electrodes in our simulations. 

The edge carbon atoms will catch any bond partner 
they can possibly get to saturate their dangling bonds. 
We assume that ZGNR edges are passivated by hydrogen, 
while edge atoms of the nanopore can be bonded cova- 
lently to either hydrogen (H-pore) or nitrogen (N-pore). 
Prior to transport calculations, we use DFT to relax the 
coordinates of all atoms within finite-ZGNR+nanopore 
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FIG. 4: The zero-bias electronic transmission function Eq. 
for an infinite homogeneous 14-ZGNR, whose edge carbon 
atoms are passivated by hydrogen, and the same nanorib- 
bon with empty H-pore or N-pore of diameter ^1.2 nm (see 
Fig. [T} drilled in its interior. 




FIG. 5: The self- consistent Hartree potential at zero bias volt- 
age for the central region of 14-ZGNR biosensors (Fig. [I]) with: 
(a) empty H-pore; (b) H-pore with cytosine positioned in its 
center within the T/z-plane [Fig. |3je)]; (c) empty N-pore; and 
(d) N-pore with thymine positioned in its center within the 
yz-p\<me [Fig.^e)]. 



or finite-ZGNR+nanopore+nucleotide until the forces on 
individual atoms are minimized to be smaller than 0.05 
eV/A 2 . The converged result of this procedure is illus- 
trated in Fig. |3jb)-(d) which shows how carbon and hy- 
drogen atoms around the nanopore move away from it so 
that the edge of ZGNR acquires a slight curvature. 

The early theoretical studies of ZGNR-based de- 
vices have utilized a simplistic tight-binding model [26] 
with single ir orbital per site and nearest neigh- 
bor hopping only, or its long-wavelength (continuum) 
approximation — the Dirac-Weyl Hamiltonian [27] — valid 
close to CNP. However, to make connection to real- 
istic device applications requires to take into account 
charge transfer [24] between different atoms [28] that 
can be used to passivate edges or chemically function- 
alize graphene, as well as the charge redistribution [24] 
when finite bias voltage is applied. For example, the 
tight-binding model with the nearest-neighbor hopping 
predicts [XT] [26] that zero-temperature conductance of 
an infinite homogeneous ZGNR is G = 2e 2 /h around the 
CNP and that current density profile is peaked in the 
middle of ZGNR despite transverse part of the eigen- 
functions having maximum around the edges [T7] . 

On the other hand, first-principles methods find that 
the zero-temperature conductance of an infinite homo- 
geneous ZGNR is G = 6e 2 /h around the CNP while lo- 
cal current is confined to flow mostly around the zigzag 
edges [20 . This is illustrated by quantized steps in the 
transmission function in Fig. [4] where T = 3 around the 
Fermi energy E — Ep = 0, and the zero-temperature 
conductance is given by the simplified version of Eq. pi), 

9 2 — 

G = \T{E). In the absence of any DNA base, the 
transmission function T(E) plotted in Fig. [4] remains 



large T — 2 around CNP E — Ep = for an infinite 
ZGNR with either H-pore or N-pore. This finding cor- 
roborates our conjecture that nanopore in the interior 
of a ZGNR is not able to substantially modify the cur- 
rent flow inherited from a homogeneous nanoribbon since 
the local current density is mostly confined around the 
edges for electrons injected at energies sufficiently close 
to E — E F = 0. We note that using spin-unrestricted 
DFT reveals the presence of edge magnetic ordering and 
the corresponding band gap opening in ZGNRs which, 
however, is easily destroyed at room temperature [29] [30] 
so that for realistic device operation ZGNRs can be con- 
sidered to be metallic [30] . 

The change in the room-temperature conductance 
of empty nanopores in Fig.[2|a) and nanopores with 
inserted nucleotide in Fig. ^Fb) is more pronounced 
when the pore is terminated with nitrogen. Since re- 
liability of predictions of NEGF-DFT simulations re- 
quires careful selection of the basis set and pseudopo- 
tentials in the DFT part of the calculation [31], Fig. [2] 
plots conductances obtained using two different compu- 
tational implementations of the NEGF-DFT formalism. 
Our home-grown MT-NEGF-DFT code [15] utilizes 
ultrasoft pseudopotentials and Perdew-Burke-Ernzerhof 
(PBE) parametrization of the generalized gradient ap- 
proximation (GGA) for exchange-correlation functional 
of DFT. The localized basis set is constructed from atom- 
centered orbitals (six per C atom, four per H atom, 8 per 
N atom, and 8 per O atom) that are optimized variation- 
ally for the electrodes and the central molecule separately 
while their electronic structure is obtained concurrently. 
For comparison, we also used commercial ATK code [16] 
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where pseudoatomic local orbitals are single-zeta polar- 
ized on C and H atoms and double-zeta polarized on N 
and O atoms. In the case of ATK, we use Troullier- 
Martins norm-conserving pseudopotentials and Perdew- 
Zunger (PZ) parametrization of the local density approx- 
imation (LDA) for the exchange-correlation functional of 
DFT. Importantly, both first-principles simulations yield 
very similar results for the conductance, as demonstrated 
in Fig. [2] 

To explain the mechanisms by which DNA bases mod- 
ulate charge transport in a ZGNR with a nanopore, we 
plot in Fig. [5] the self-consistent Hartree potential within 
the central region of our biosensor at zero bias voltage 
obtained by solving the Poisson equation with the bound- 
ary conditions that match the electrostatic potentials of 
two attached ZGNR electrodes. We see that there is 
a substantial difference in this potential when switching 
from an empty pore to nanopore containing a DNA nu- 
cleotide. In the examples in Fig. [5| cytosine is inserted 
into the H-pore and thymine into the N-pore — these are 
the situations for which there is the largest change in con- 
ductance in Fig. [2] when compared to the corresponding 
empty nanopores. 

An important issue [TJ |9] for the uniqueness of the 
conductance modulation signal associated with each nu- 
cleotide is to examine how such signal gets modified 
when varying the orientations of the nucleotide within a 
nanopore. For selected orientations shown in Fig. [3^b)— 
(e), the conductance variation for all four DNA nu- 
cleotides is plotted in Fig. [3ja). We find small overlap 
between conductance distribution for T and A or A and 
C, and no overlap between conductance intervals for T 
and C or C and G. Nevertheless, the intervals in Fig.|3ja) 
should be considered as setting only the limits on con- 
ductance variation since not all values within the interval 
will be sampled experimentally. That is, some of the nu- 
cleotide positions in Fig.[3^b)-(e) are selected to generate 
maximum conductance variation, and they would require 
significant bending of the DNA molecule to put the nu- 
cleotide into such position with respect to the nanopore. 
Instead, to find the most probable fluctuations in orien- 
tations of DNA nucleotides, the NEGF-DFT calculations 
of conductance should be coupled to molecular dynamics 
simulations of DNA translocation through the nanopore 
including hydrodynamic interactions with the surround- 
ing solvent [32] . 

Finally, in Fig. [6] we clarify the range of operating bias 
voltages that ensures a linear-response regime for our 
biosensor where the measured current is given simply by 
multiplying conductances in Figs. [2] and [3] by the bias 
voltage. Both current-voltage characteristics in Fig. |6j 
computed for a biosensor with an empty N-pore and the 
same pore containing guanine, behave linearly within the 
interval ~ -0.05 V to ~ 0.05 V. 

In conclusion, using first-principles quantum transport 
simulations, we investigated a novel type of graphene 
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FIG. 6: Current- volt age characteristics of 14-ZGNR with N- 
pore which is empty (dashed line) or contains guanine (solid 
line) in its center placed within the yz-p\&ne [Fig. |3je)]. The 
current at finite bias voltage is computed using our home- 
grown MT-NEGF-DFT code [T31[T5] . 



nanopore-based sensors for rapid DNA sequencing which 
rely on nucleotide-specific modulation of a large trans- 
verse conduction current (of the order of fiA at bias 
voltage ~ 0.1 V). This is achieved by exploiting unique 
features of the electronic transport through graphene 
nanoribbons with zigzag edges where local current den- 
sity is confined mostly around the nanoribbon edges. An- 
other candidate nanowire carrying edge currents are re- 
cently fabricated [18 chiral GNRs. Thus, the nanopore 
in the GNR interior cannot substantially diminish the 
edge currents, whose magnitude is then modulated by 
the passage of nucleotides in the course of DNA translo- 
cation through the pore. Our analysis demonstrates that 
each DNA nucleotide will generate a unique electrostatic 
potential that modulates the charge density in the sur- 
rounding area. The operating current, which is several 
orders of magnitude greater than the tunneling current 
employed in previously considered biosensors with trans- 
verse electron transport [6HT3] is expected to be much 
larger than its fluctuations due to thermal vibrations of 
the graphene membrane or structural fluctuations of the 
translocated DNA molecule. Furthermore, the device 
remains in the linear-response regime for bias voltages 
< 0.05 V. 

We thank S. Sanvito for illuminating discussions. This 
work was supported by DOE Grant No. DE-FG02- 
07ER46374 (K. K. S. and B. K. N.) and NIH Grant No. 
R21HG004767 (M. D.). The supercomputing time was 
provided in part by the NSF through TeraGrid resource 
TACC Ranger under Grant No. TG-DMR100002. 



6 



* Electronic address: bnikolic@udel.edu 

[1] E. E. Schadt, S. Turner, and A. Kasarskis, Human Molec- 
ular Genetics 19, R227 (2010). 

[2] C. A. Merchant et al, Nano Lett. 10, 2915 (2010). 

[3] G. F. Schneider et a/., Nano Lett. 10, 3163 (2010). 

[4] S. Garaj et al, Nature 467, 190 (2010). 

[5] A. K. Geim, Science 324, 1530 (2009). 

[6] H. W. C. Postma, Nano Lett. 10, 420 (2010). 

[7] J. Prasongkit et a/., Nano Lett. 11, 1941 (2011). 

[8] Y. He et a/., Adv. Funct. Mater. 21, 2674 (2011). 

[9] T. Nelson, B. Zhang, and O. V. Prezhdo, Nano Lett. 10, 
3237 (2010). 

[10] M. Zwolak and M. Di Ventra, Rev. Mod. Phys. 80, 141 
(2008). 

[11] V. Meunier and P. S. Krstic, J. Chem. Phys. 128, 041103 
(2008). 

[12] M. Tsutsui, M. Taniguchi, K. Yokota, and T. Kawai, 

Nature Nanotech. 5, 286 (2010). 
[13] S. Huang et a/., Nature Nanotech. 5, 868 (2010). 
[14] K. K. Saha et a/., Phys. Rev. Lett. 105, 236803 (2010). 
[15] K. K. Saha, W. Lu, J. Bernholc, and V. Meunier, J. 

Chem. Phys. 131, 164105 (20 09). 
[16] http://www.quantumwise.com 

[17] L. P. Zarbo and B. K. Nikolic, Europhys. Lett. 80, 47001 



(2007). 

[18] C. Tao et a/., Nature Physics 7, 616 (2011). 
[19] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 

[20] D. Areshkin and C. White, Nano Lett. 7, 3253 (2007). 

[21] J. Cai et a/., Nature 466, 470 (2010). 

[22] X. Jia et a/., Science 323, 1701 (2009). 

[23] J. Taylor, H. Guo, and J. Wang, Phys. Rev. B 63, 245407 
(2001); M. Brandbyge, J.-L. Mozos, P. Ordejon, J. Tay- 
lor, and K. Stokbro, Phys. Rev. B 65, 165401 (2002). 

[24] D. A. Areshkin and B. K. Nikolic, Phys. Rev. B 81, 
155450 (2010). 

[25] Y. He et al, Appl. Phys. Lett. 97, 043701 (2010). 

[26] A. Rycerz, J. Tworzydlo, and C. W. J. Beenakker, Nature 
Phys. 3, 172 (2007). 

[27] L. Brey and H. A. Fertig, Phys. Rev. B 73, 235411 (2006). 

[28] F. Cervantes-Sodi, G. Csanyi, S. Piscanec, and A. C. 
Ferrari, Phys. Rev. B 77, 165427 (2008). 

[29] O. V. Yazyev and M. I. Katsnelson, Phys. Rev. Lett. 100, 
047209 (2008). 

[30] J. Kunstmann, C. Ozdogan, A. Quandt, and H. Fehske, 

Phys. Rev. B 83, 045414 (2011). 
[31] M. Strange, I. S. Kristensen, K. S. Thygesen, and K. W. 

Jacobsen, J. Chem. Phys. 128, 114714 (2008). 
[32] M. Fyta, S. Melchionna, and S. Succi, J. Polymer Science 

Part B: Polymer Physics 49, 985 (2011). 



